%% Residency times within the corrugated and straight channels
clear all; close all;
%%%%%%%%%%%%%%%%%%%% root
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
root = 'D:\';
%%%%%%%%%%
load([root, 'Compilation DiodeNatCom\figure 3\residency.mat']);

hFig = figure('color',[1,1,1],'Units','centimeters','Position',[5,5,1.2*4.5,1.2*4]); hold on; box on;
cmap =[[0,0,1];[1,0,0];[0,0.75,0];[140,100,0]/255];

for r=1:4
plot(residency(r).x_L,residency(r).y_L,'color',cmap(r,:),'LineWidth',1);
plot(residency(r).x_R,residency(r).y_R,'color',cmap(r,:),'LineWidth',1);
end

yticks([1,10,100,500]);
ylim([0.3,500]);
xlim([0.05,10]);

set(gca,'YScale','log');
set(gca,'XScale','log');

xlim([0.1,10])

